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Abstract 

Lambda flow in Ni+Ni collisions at SIS energies is studied in the relativistic 
transport model (RVUU 1.0). It is found that for primordial lambdas the 
flow is considerably weaker than proton flow. Lhe inclusion of final-state 
interactions, especially the propagation of lambdas in mean-field potential, 
brings the lambda flow close to that of protons. An accurate determination of 
lambda flow in heavy-ion experiments is shown to be very useful for studying 
lambda properties in dense matter. 
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I. INTRODUCTION 



Collective flows of hadrons and light fragments in heavy-ion collisions have been observed 
unambiguously at various incident energies JT| — ^] . Detailed comparisons of theoretical pre- 
dictions with the experimental data have already provided valuable information about the 
nuclear equation of state IM, the in- medium nucleon-nucleon cross sections [TD|, and the 



medium modification of hadron properties [IT]. It has been found that proton and fragment 
flows are largely caused by the compressional pressure generated at high densities and thus 
carry the information about the nuclear equation of state. Proton flow at intermediate en- 
ergies has also been found to be sensitive to medium modifications of the nucleon-nucleon 
cross sections. Transport model calculations without including mean-filed potentials have 
shown that pion flow fll2"l-|l4|j , as well as antikaon and antiproton flow [IJ||n|, are mostly 
determined by their large absorption cross sections by nucleons, and thus usually show 'an- 
tiflow' with respect to nucleons. These results are, however, qualitatively changed when 
their mean-field potentials are included in the transport model. Both the antiflows of pions 
17j and antikaons |18[ are found to be significantly reduced by their attractive potentials 



in nuclear medium. A similar effect is expected for the antiflow of antiprotons when their 
attractive mean-field potential is included in the study. Also, kaon flow has been shown 



to be sensitive to the kaon potential in nuclear medium 1TTJ. A comparison with recent 
experimental data from the FOPI collaboration at GSI has shown that a kaon feels both 
an attractive scalar and a repulsive vector potential [|l9j, consistent with the prediction of 
effective chiral Lagrangians [pOj^T| . The resulting kaon potential is weakly repulsive and 
thus makes the kaon flow extremely weak compared to nucleon flow. 

More recently, collective flow of lambdas has been measured through the reconstruction 
of pn~ pairs by the FOPI collaboration in Ni+Ni collisions at 1.93 GeV/nucleon ||, and by 
the EOS collaboration in Ni+Cu collisions at 2.0 GeV/nucleon Although the magnitude 
of lambda flow differs in the two measurements, both collaborations have found that lambdas 
flow in the same direction as nucleons, and its strength is similar to that of proton flow. The 
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quantitative difference between the two measurements is due to different acceptance cuts 
and centrality selections in the two experiments. 

In heavy-ion collisions around 1-2 GeV/nucleon, lambdas are mainly produced in as- 
sociation with kaons from baryon-baryon and meson-baryon collisions. In Refs. [0,0 we 
have studied kaon flow from Ni+Ni collisions at similar energies, and found that the flow of 
primordial kaons is considerably weaker than proton flow. This is also true for primordial 
lambdas, as we shall show later. In the case of kaons, the inclusion of final-state interactions, 
mainly their propagation in mean field, tends to repel them away from nucleons, and this 
reduces kaon flow or even changes its sign if the kaon potential is very repulsive. 

The final-state interactions of a lambda are different from that of a kaon. At low energies, 
the lambda-nucleon cross section is relatively large compared to the kaon-nucleon cross 
section. Also, instead of a weak repulsive potential as for a kaon, the lambda potential 
is attractive in the density and momentum range relevant for heavy-ion collisions at beam 
energies around 1-2 GeV/nucleon. These differences are expected to make the final lambda 
flow quite different from the kaon flow, although they are very similar without final-state 
interactions. In particular, we shall show that the attractive potential between a lambda 
and the nuclear matter makes its flow similar to the nucleon flow. 

The purpose of this paper is to study quantitatively the effects of final-state interactions 
on lambda flow in heavy-ion collisions at SIS energies. We also study the sensitivity of 
lambda flow to its properties in dense matter, which are important for understanding the 



properties of 'strange' stars and the possible kaon condensation ]22|- g4j] . In Section II, we 
briefly review the relativistic transport model (RVUU 1.0) that has been used extensively 
in studying heavy-ion collisions at SIS energies. We also discuss how we treat lambda 
production, scattering, and propagation in the transport model. The results from this study 
are presented in Section III. The paper ends with a short summary in Section IV. 
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II. THE RELATIVISTIC TRANSPORT MODEL 



The relativistic transport model (RVUU 1.0), first developed in Ref. p5 and further 



extended in Refs. [|6|,^7], is based on the non- linear o-uj model and includes explicitly 
the nucleon, delta resonance, and pion. It can also treat the production of eta, kaon, 
antikaon, hyperon, antiproton, and dilepton in heavy-ion collisions at SIS energies using the 
perturbative test particle method. 

At SIS energies, lambdas are mainly produced from baryon-baryon [BB — ► BAK, with 
B denotes either a nucleon or a delta resonance) and pion-nucleon (irN — > AK) collisions. 
The process 7rA — > AK is neglected since its cross section is much smaller than that of 
ttN — > AK, as shown in the resonance model of Ref. p8[. For the isospin-averaged cross 



sections (Jbb^bak, we use the parameterization of Ref. [p9[1 , e.g., 



with 

Pr, 



2v^ 



(Tnn-,nak = 0.072^^ mb, (1) 
m K 



1 l/ 2 

(s - (m N + m A + m K ) 2 )(s - (m N + m A - m K ) . (2) 



The cross section for ttN — > AK is taken from the parameterization in Ref. pO 



i.e., 



a*N->AK = 4.94(v^ - yfa) mb, - v 7 ^ < °- 091 GeV ( 3 ) 
U ' U J ° mb, v^- 0.091 GeV, (4) 



0.01 + (v^~ yfid 
with y/so = m\ + m K . 

In the following calculations, we always include the medium modification of kaon proper- 
ties by using its in-medium (pole) mass m* K determined from the mean-field approximation 
to the chiral Lagrangian f26|, i.e., 



* n ^KN 3 1/2 

m K ^m K [l-—^p s + — p N 7 , 5 

f 2 m z K 4:f 2 m K 

where / = 93 MeV is the pion decay constant, and Y^kn ~ 350 — 450 MeV is the KN 
sigma term, which depends on the strangeness content of a nucleon. The scalar and nuclear 
densities are denoted, respectively, by p s and p^. 



The lambda mean-field potential in nuclear medium has been extracted from the prop- 



erties of hypernuclei [|31],0] . It has also been studied in the Dirac-Brueckner approach using 
a boson-exchange model for the AiV potential [Q, or by extending the Walecka model from 
SU(2) to SU(3) [0,[3^]. In the naive constituent quark model, the lambda mean- field po- 
tential is about 2/3 of that of a nucleon, and is in qualitative agreement with the empirical 
observation |H],0]. To see the sensitivity of lambda flow to the lambda potential in dense 
matter, we adopt the constituent quark relation but introduce a parameter a in the lambda 
vector potential, i.e., 

E£«2/3E£, S#«2/3aX#, (6) 
where E^ and Ey are the nucleon scalar and vector potentials in the non-linear a-u model 



|27|. The lambda optical-model potential can then be defined as 

U A (p, p) = ((m A - E* ) 2 + p 2 ) 1/2 + E* - ( m \ + ^ . (7) 



We consider three values of a, i.e., 0.85, 1.0, and 1.2 so that the lambda potential is varied 
within a reasonable range. The corresponding lambda potential is shown in Fig. 1. The 
nucleon potential used in this study is based on the parameter set in Ref. [0] that corre- 
sponds to a soft equation of state with a compression modulus of 200 MeV and a nucleon 
effective mass of 0.83mjy at normal nuclear matter density. The solid circle with error bars 
is the currently determined lambda potential in nuclear matter from both the structure of 
hypernuclei [0,0] and the Dirac-Brueckner calculation |33| . 



Including the medium modification of the lambda, as well as the nucleon and kaon, y/s in 
Eq. (D is replaced by \/s* = (m^ + p 2 ) 1 / 2 + (m$ 2 + p 2 ) 1 / 2 + (1 - 2/3a)E^, while m N , m A , 
and rriK are replaced by m* N = — Tig , m\ = — Eg, and m* K , respectively. Similarly, 
y/s and y/s^ in Eq. (|) are replaced by y/s* = (m% + p 2 ) 1/2 + (m 2 + p 2 ) 1 / 2 + (1 - 2/3a)E^ 
and a/so = rn* A + m* K , respectively. In calculating lambda production, we have thus included 
the change of the threshold due to the medium modification of hadron properties. 

After production, a lambda interacts with surrounding baryons. This includes both AA^ 
scattering and lambda propagation in the mean field generated by nuclear medium. In 



principle, the AA" scattering is modified in medium, but we have not attempted to include 
this effect in a consistent fashion. As in most transport models, the free cross section will 
be used. However, we will show results from varying the AA" cross section. Furthermore, 
we will not address the question of consistently separating the interactions into a mean-field 
and a scattering part in the transport model simulation, which has recently been raised in 



Ref. [36 



The AA" cross section has been measured experimentally |37j and also studied in a boson- 
exchange model ||33|| . For beam energies considered in the present study, the AA" inelastic 
scattering, mainly A A" — > A A, is relatively unimportant and is neglected. We thus consider 
only AA" elastic scattering. The experimental data can be fitted by 

cr AN = 12.0 + ^ mb, (8) 

Plab 

where pi a b, in units of GeV, is the momentum of the lambda in nucleon rest frame. This cross 



section is shown in Fig. 2, where circles are data from Ref. |H] and the solid curve is the 
parameterization. The angular distribution for the AA" scattering is assumed to be isotropic 
in its center-of-mass frame, which is consistent with the results from the boson-exchange 
model of Ref. . To see how the lambda flow might depend on the angular distribution in 
AA" scattering, we have carried out a calculation by taking the AA" differential cross section 



as the proton-proton one |p8| , and we find that this has very little effect on the final lambda 
flow. 

Including a mean-field potential, the lambda equations of motion between collisions are 
then given by 

dx p* dp „ rT , 



where E* = Jp* 2 + (m A — S^) 2 . 
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III. RESULTS AND DISCUSSIONS 



Using the model described in the above, we have carried out a calculation of lambda 
production in Ni+Ni collisions at 1.93 GeV/nucleon and impact parameter b < 4 fm that 
corresponds approximately to the centrality selection in the FOPI experiment. A direct 
comparison of the predicted lambda flow with the preliminary FOPI and EOS data requires 
also knowledge on the experimental acceptance cuts. Since these corrections apply similarly 
to protons and lambdas, we have not attempted to include them in this study. What is 
significant in both FOPI and EOS data is the fact that the observed lambda flow is in the 
same direction as proton flow, and both have similar strengths. We thus compare instead 
the predicted lambda and proton flows from the transport model. 

The results for lambda flow, i.e., the average transverse momentum (p x ) as a function 
of rapidity y in the nucleus- nucleus center-of-mass frame, are shown in Fig. 3. The dashed 
curve is the flow of primordial lambdas. The dotted-dashed curve gives the flow of lambdas 
after including elastic AN scattering with a cross section given by Eq. The lambda flow 
including both scattering and lambda propagation in mean field (with a = 1.0) is shown by 
the solid curve. For comparison, we also show by dotted curve the proton flow. 

To be more quantitative, we introduce the flow parameter defined as the slope of the 
average transverse momentum at mid-rapidity, i.e., 

v d(Px) 



\v=o> 



(10) 



dy 

As in the case of kaons, the flow of primordial lambdas, which arises mainly from Lorentz 
boost in the direction of the baryon-baryon or meson-baryon pairs that produce the antikaon, 
is considerably weaker than the proton flow. Their flow parameter is about 60 MeV as 
compared to the proton flow parameter of about 140 MeV. Including elastic AN scattering 
increases the flow parameter to about 80 MeV. This is due to the effect of thermalization that 
increases the magnitude of lambda momentum. The propagation of lambdas in mean-field 
potential further enhances their flow in the direction of nucleons, and the flow parameter is 
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now about 115 MeV, and is close to that of protons. Overall, final-state interactions enhance 
the lambda flow parameter by about a factor of two. 

The importance of final-state interactions can also be seen from (p x ) ma x, which occurs 
near the projectile and target rapidities. The magnitude of (p x )max for primordial lambdas 
is about 45 MeV, which increases to about 60 MeV when AiV scattering is included, and 
further increases to about 100 MeV after including also lambda propagation. Final-state 
interactions thus increase {p x )max by about a factor of two as well. 

To learn quantitatively the lambda potential from lambda flow requires not only more 
accurate experimental data but also better theoretical understandings on how the flow is 
affected by a possible change of AiV cross section in nuclear medium. For this purpose, we 
have carried out two calculations, one with a aat reduced by 50% and the other with <j\ N 
increased by 50% (a is always kept at 1.0). The results are shown in Fig. 4. It is seen that 
changing AA" cross section by a factor of 3 does not affect significantly the lambda flow, with 
the flow parameter modified by only about 15 MeV. On the average, each lambda undergoes 
about 1.3, 2.5 and 3.6 collisions for 0.5<tatv, cr\ N , and 1.5a\N, respectively. This implies that 
with 1.3 AA" scattering, lambdas are already thermalized by nucleons, so further increase of 
the cross section does not have significant effects on the lambda momentum distribution. 

To see the sensitivity of lambda flow to the lambda potential we have done two additional 
calculations with a = 1.2 and 0.85, while always using the AA" cross section given by Eq. 

The results for lambda flow are shown in Fig. 5 together with that using a = 1.0. We 
see that with an increasingly attractive lambda potential, the lambda flow becomes stronger. 
The flow parameter F changes from about 95 MeV for a— 1.2 to about 135 MeV for a=0.85. 
It is thus possible to differentiate the scenario with a shallow lambda potential (as given by 
a=1.2) from that with a deep lambda potential (as given by a=0.85) if the lambda flow can 
be measured with good accuracy. 

Similarly, one can also look at the azimuthal distribution of lambdas to examine the 
effects of lambda potential. The results for the lambda azimuthal distribution near the 
target rapidity are shown in Fig. 6 for three different potentials, as well as for the case 
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without lambda potential. We have normalized these results around <p = 0°. In all four 
cases, the lambda azimuthal distribution exhibits a peak around = 180°, as in the nucleon 
distribution shown in the figure by the dotted curve. The lambda azimuthal anisotropy 
gets more pronounced as the lambda potential becomes more attractive. To show this more 
quantitatively, we introduce an anisotropy parameter R defined by 

dN '/ 'dxf>{<f> = 180° ± 15°) 

dN/d<j>(<j) = 0° ± 15°) ' (U) 

It is seen that R changes from about 1.9 to about 2.6 when a changes from 1.2 to 0.85. 
Without lambda potential R is about 1.4. For comparison, we note that the anisotropy 
parameter for nucleons is R ~ 3.0. Thus, the lambda mean-field potential also enhances its 
azimuthal anisotropy in the direction of nucleons. 

In Fig. 7, we summarize the dependence of the flow parameter F and the anisotropy 
parameter R, as well as the lambda yield, on the strength parameter a of the lambda vector 
potential. It is seen that the lambda yield also depends on the strength of lambda potential. 
A stronger repulsive lambda vector potential increases the threshold for its production and 
thus decreases its yield. As a result, we see a similar dependence of the lambda yield on the 
strength parameter a as for the lambda flow parameter and anisotropy anisotropy. 

IV. SUMMARY 

In summary, we have studied lambda flow in heavy-ion collisions at SIS energies. We 
have found that the primordial lambdas show a relatively weak flow as compared with 
the nucleon flow. The inclusion of final-state interactions, especially the propagation in 
mean-field potential, enhances the lambda flow in the direction of nucleons, and brings 
the theoretical results in agreement with the preliminary data from both FOPI || and 
EOS 0] collaboration. Significant differences in both lambda in-plane and out-of-plane 
flows are found between the results with and without lambda potential. On the other 
hand, the final lambda flow is relatively insensitive to changes of the AN cross section 
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within a reasonable range. Accurate measurements of both lambda yield and flow in heavy- 
ion collisions thus allow us to determine the lambda potential in the dense matter formed 
in heavy-ion collisions. This information cannot be obtained from studies of hypernuclei 
which provide only the lambda potential at and below normal nuclear matter density. For 
understanding the properties of neutron stars, lambda properties at higher densities than 
normal nuclear matter are needed. The information one derives from lambda flow in heavy- 



ion collisions is thus very useful for studying neutron star properties 22-24 
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Figure Captions 

Fig. 1: The lambda optical potential. The solid circle is the empirically determined 
lambda potential at normal nuclear matter density. 

Fig. 2: The elastic AiV cross section. Circles are experimental data from Ref. |3"7|| , and 
the curve is the parameterization given by Eq. (|5p. 

Fig. 3: The average transverse momentum of lambdas as a function of rapidity. The 
dashed, dotted-dashed, and solid curves are for primordial lambdas, lambdas with scattering, 
and lambdas with both scattering and propagation, respectively. The dotted curve is for 
protons. 

Fig. 4: Effects of AiV cross section on lambda average transverse momentum distribu- 
tion in rapidity space. The dotted, solid, and dashed curves are obtained with 0.5o"aat, oajv, 
and 1.5<7aat, respectively. 

Fig. 5: Effects of lambda potential on its average transverse momentum distribution in 
rapidity space. The dotted, solid, and dashed curves are obtained with a=0.85, 1.0 and 1.2, 
respectively. 

Fig. 6: The azimuthal distribution of lambdas near target rapidity. The dotted curve 
gives that of nucleons. 

Fig. 7: The lambda yield, flow parameter, and azimuthal anisotropy parameter as 
functions of the strength parameter a for the lambda vector potential. 
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